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Abstract. We consider in 2D the following special case of the 
Mumford-Shah functional 



J{uX} = 




It is known that if the minimizer has a cracktip in the ball Bi 
(assume at the origin), then u « AS-^/z at this point. We calculate 
higher order terms in the asymptotic expansion, where the homo- 
geneity orders of those terms appear to be solutions to a certain 
trigonometric relation. 



1. Introduction 

1.1. The Problem. Let us consider the foUowing minimization prob- 
lem; Minimize the functional 

(1) J{uJ)= ! \Vu\''dx + X\/^H\T) 

JBi\r V 2 

over all curves T G Bi starting at the point —1 of the complex plane 
and functions u G H^{Bi\T) with boundary values g on dBi. We 
will assume that g is a. little perturbation of the function X'^^/z = 
^^1/2 gin on the dBi, where the discontinuity is taken at the point 
-1. 

Let us agree that throughout the paper, depending on the con- 
text, the functions Q^/z, r" sin(a0) and r" cos{a(f)) will be defined in 
M^\{(a;, 0)| — oo < x < 0}, in -Bi\r or in Bi\rr, where the discontinuity 
is at {{x, 0)1 — oo < a; < 0}, at F or at respectively. 

It is proven in |BDj that the pair 

{\r^/^sm(f)/2]{{x,0) \ - oo < a; < 0}) 

is the global minimizer of the functional ([T]) in the plane M^. From this 
result it follows that for the boundary value function g = A sin 0/2, the 
interval Fq = {(x, 0)| — 1 < x < 0} and the function mq = Ar^/^ sin 0/2 
give the absolute minimizer. 
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We will assume the stability of the problem under small perturba- 
tions, i.e. that for small perturbations of g the minimizing set F is a 
curve connecting the point (—1,0) G dBi fl M with an unknown point 
inside the ball Bi := -Bi(O) and that this curve is not "far from" the 
interval Fq. 

The authors refer to [D] and |AFPj for further references and as 
sources of inexhaustible knowledge in the field. 

Probably some of the assumptions made in the paper could be avoided 
but still we believe that the results even in this simplified formulation 
are new and important. 

1.2. The Main Theorem. Assume the minimizer of the functional 
(|l| in 5i C is given by the pair (m,F), where F = {(-t,/(t))|0 < 
t <1}, f e C\[0, 1)) /(O) = /'(O) = 0. Further assume there exist a 
limit inC/J(0,l)) 

9it) = hm ^ 0. 

P^o maxo<r<p |j(r)| 

Then there exists a constant C 7^ such that 
f(t) = + o(t^/2+"^) 

and 

(2) u{x, y) = ASV^ + y) + CA6fcr"'= cos(«fc0) + o(r"'=), 
where k < ak < k + \ is one of the positive solutions of 

/7r a 



tanivra) 



2a^ , 



1 ' 



hk are some other absolute constants {k = 1,2, ... ), Si = and for 
k > 2 there are some constants Cj depending on u such that 



/ X sr^ . / 2j + 1 

T.k{x,y) = y cjr 2 sm 



If 



lim , ..=0 



P^o maxo<r<p 
then limt^o t-'^^f{t) = for any M > 0. 

1.3. Acknowledgment. The authors are grateful to MSRI at UC 
Berkeley for the invitation to attend the Free Boundary Problems Pro- 
gram in 2011 and particularly to Henrik Shahgholian, who was the 
speaker of the program. The substantial part of the research was car- 
ried out during the stay in Berkeley. The second author thanks Stephan 
Luckhaus for useful discussions. 
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2. Asymptotic of the minimizing function near the 

cracktip 

In this section we want to calculate the asymptotic expansion for the 
minimizer function u near the cracktip. 

Let us assume that pair {u, T) is as above and minimizes the func- 
tional ([T]). Without loss of generality we can assume that the origin 
is the endpoint of F and is tangential to the x-axis at that point, for 
this reason we allow F to "start" from a point different from (—1,0). 
We will parametrize F = {{—t,f{t)),t G [0, 1]} and at the cracktip we 
have /(O) = /'(O) = 0. 




We know that on F the balance condition is fulfilled (see |AFPj ). i.e. 

(3) x'Hr{~t,m) = ^'jYTjnW ^ 

where [|Vm(2;)P]^ = \Vu{x,y + 0)1"^ — |Vm(x, y — 0)p is the difference 
of the values of \Vu\^ on the two sides of the discontinuity set and Ht 
is the curvature of F. Further we know that the first asymptotic term 
of the function u near the origin is Ar^/^ sin 0/2. A natural question 
arise what are the next terms in the asymptotic? 
Let us consider the function 

(4) vp{x,y) = S;\uipx,py)-p'/'\r'/'sm{(P/2), 

where Sp = supBp\u{x,y) — Ar^/^ sin(0/2)| and Spp~^^'^ — )■ as p — )■ 0. 
Note that 

r^/^sin(0/2) = -^sgnp(x, y)\/ ^/x^ + y^ ~ x, 

where sgnp(x, = sgn(?/) if x > 0, sgn-p{x,y) = 1 if x < 0,y > /(— x) 
and sgnp(x,?/) = —1 if x < 0,y < /(— x). 
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Thus 



VrV'sin((/./2) = I -1 + 



X 



2V2sgnj,{x, y)\y\ \ ^/x^ + ^2' ^x'^ + y'^^ 

The rescaled crack Tp is given by {(— fp{t))} where fp{t) — p~^f{pt) 
Let us calculate 9,^fp and drVp on the crack, where the vectors r = 

(1 + f'lrH-l; f'p) and i/ = (1 + flrHfp; 1) are the tangential and 
normal directions on Fp 

Since di,u = on F we have that 



(5) d,Vpi-t, fp{t)) = V ■ Vr'/^ sin(0/2) = 



/' -1 - 



+ 



/p 



I 



f 



\ 



-1 - 



v 



J 



If we now denote by ap = supte{o,i)\fpit)\ = p ^supte(o,p)\f{t)\ and 
Qp — ^we obtain that 



(6) d,vp{-t,fp{t)) 



P 



SAVi 



( 



\ 



\ 



1 + 



9p 



1 + 



In the Appendix A at the end of the paper we prove that liminfp^o J? < 

p ^ Sp 

+00. Assume for a sequence — >■ 

(7) ^ ^ > 

Pk ^ ^pk 

and pass to the limit taking if necessary a subsequence of pk to have 
9pk ~^ 9o and Vp^, vq] we obtain that for x < 
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dyVo{x,±0) = =F- 



XA 



-X 



9o 



-X 



Let us now rewrite the condition ([3]) in terms of / and Vr- After 
rescaling we have 



(9) 



2i± 



The mean curvature can be expressed in terms of function / 



H{T,){-tJ,{t)) 



(1 + im''^ 



and Ot-u^ in terms of Vp 



r ■ Vr^/^ sin(0/2) can be expressed in ters of / as well 



(10) [r-V//^in(0/2)](-t,/,(t)) 



(1 + /1)-A-^' 



t 



1 + 



vv 



(l + o(/('Vt)) 



fp 



+ 



fpf p 



fp 



4v^ 



VV 



fp 



fpf 



\ 



pj p 



J 



fp 



J 



(i + o(/Vt))^(i + o(/pA)) 
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Putting together we rewrite ^ as follows 



(11) X'pa, 



(1 + im'/' 



^[{SAM-t, m) ± //'A [r • sin(0/2)](-t, /.(t))'' 



2^,Ap^/2(l+o(/(')/t))^(l+o(/,/t))(9.<(-t,/,(t)+9.t;;(-t,/,(t) 



Thus 



(12) A^—""^ 



(1 + (/;)2)3/2 

2A(1 + o(/(')A))^(i + o{f,/t)) {d^v^i-t, f,{t) + d^v;{~t, f,{t)) 
and passing to a limit 



(13) AA(7;'(-x) = {-x)-'/'^-[d,vUx,0) + d,v^{x,0) 
Let us now denote by 



Wi{x,y) := -(t;o(x,y) - Vo{x, -y)) 



and by 

Then we obtain 



W2{x, y) := -{vo{x, y) + Vq{x, -y)). 



2dyWi{x,Q) = 0, for X < 0, 



which means Wi {x, y) is a harmonic polynomial of the form '^z 

ken. 

For W2 we obtain 
(14) 2dyW2{x, 0) = {-2g',{-x) + '^"^""^ 



2fc+l 
2 , 



-X 



and 



2(-a:)-^/^y|a.iy2(x,0) = AA^7[;(-x). 

If y4 = then W2 = 0, and then we have Vo{x,y) = Wi{x,y) = 
in M2\{(x, 0)1 - 00 < X < 0}, A; e N. We will come back to this 
case at the end of this section. 
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Let us now consider the more interesting case A > 0. 

First asume W2 = 6r"cosa0, (we know that a > 1/2). We obtain 



9 [r) = ' cosTra 



and 



9{r) 



2b [nacosTTa 

r 2 . 



V 2 a2 _ 1 



From here we first obtain, since sup(g(o g{t) = 1 and g{r) = r""'"2 



(15) 



2b 



Tx a cos vra 
2^CT 



and then from (14) 



-2asin7rar" ^ = "^^^ f — 2 , „ , 



(a + ^)r"-^/' + r"-i/2 j = -aAAr""^ 



thus 



or 



(16) 



smyra 



tan Tra 



TT a cos vra 
2 a2_ 1 



2a2 _ 1- 



Thus we obtain that a = a^, where < < + ^ are the positive 



solutions of (16). Note that the smallest positive solution ai ~ 1.2739. 




Picture [2] 
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We proved that if we assume that 

W2{x,y) = ±r"cos(a0). 
then a = ak for some k = 1,2, . . . and 

g{t) = ±ti/2+a_ 

In the Appendix B we prove that the case W2 = br" cos a(j) is the only 
case we need to consider. 

Since the orders of possible solutions for Wi and W2 are different we 
will always have that one of them must be identically zero. This means 
that 6 = ±1 and from ([T]) we obtain that if 

f{t) = cr'^+^z^ + o(t"*+i/2^ 

then 

CX 

u{x, y) = X^^/z ± -pr'""'" cos{ak(p) + o(r"'=), 



and we see from (15) that = ^ > depends only on k 



Let us now come back to the case A = 0: that means W2 vanishes 
and Wi is not. If we have that 



u{x,y) = \'^^/z + J:k{x,y) + o{r 2 ), 

where Sfc is as in (|2]), then we can repeat the arguments above with 
the function 

Vp{x, y) = S^^{u{px, py) - p^/^Ar^/^ sin(0/2) - Sfe(pa;, py)), 

where Sp = supBp\u{x,y) — Xr^^'^ sin(0/2) — T,k{px, py)\, instead of Q. 
This means that we can iterate this procedure until we get the first 
non-zero W2 term. If the iteration goes to infinity without giving any 
non-vanshing W2, then we have that t~^^ f{t) — )■ as t — )■ for any 
M > 0. Which would correspond the case lim„ -.n s, = 0. 

f p->u maxo<r<p \fij)\ 

3. The missing Euler-Lagrange condition 

In this section we present some heuristic arguments to justify the 
existence of a missing Euler-Lagrange condition. 

As it is proven in |BDj . for the boundary function g = X sin 0/2, the 
interval F = {(x, 0)| — 1 < a; < 0} and the function u = Ar^/^ sin</)/2 
give the absolute minimizer of ([T]). Let us take a point w E Bs and 
minimize the functional ([T]) among all curves starting at the disconti- 
nuity point at the boundary data g, i.e. (—1, 0) and ending at the point 
w. Heuristically it is natural to assume that for small 6 the minimizers 
{uw, Tui) will exist and F^^, will be a smooth curve. This would mean 
in particular that the first term in the asymptotic of Um will be, up 
to a rotation, of the form Xw'^\^¥^^w. The authors think that Xy^ is 
monotone in the horizontal direction and this would mean that there 
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exists an interface S G Bs on which A^, = A, where A is the "right" 
constant from the functional ([T]). 

If we now look at the pairs {u^, F^), where w G S, then all of them 
satisfy the three known Euler-Lagrange conditions 

(i) Au^ = in Bi\T^ 

(ii) X^Hr^ = ^[\Vu^\Y on 

(iii) A^ = A 

On the other hand we see that only one of them is the absolute 
minimizer. Since S (if the heuristics is true) is an interface near the 
origin the missing condition should be of first order. 

From the Main Theorem one can see that the curvature at the 
cracktip can be either zero or oo (in case the coefficient of the ai-term 
is not zero). One possible candidate for the missing Euler-Lagrage 
condition could be the following: 

The curvature of the crack vanishes at the cracktip. 

or in terms of asymptotic, that the coefficient C of the r"^ cos(q;i<^)- 
term in the expansion (|2| vanishes. 

Appendix A 
Here we will prove that the expression 



is bounded at least on a subsequence near the origin. 

To make the problem easier to handle let us "open" the crack. We 
consider as complex plane and apply ^/z conformal transformation 
given by 

/ \ , f^S^r{x,y) / / 2 , 2 I sgnr(x,y) / / ^ , 2 ^ 
{x,y) ^ ( ^ ^/x^ + + x; ^ ^/x^ + - x j . 
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The new function w{z) := u{z'^) will be harmonic in fl, where fl 
is the new domain. The free discontinuity set F will be mapped to 
the symmetric set S = dfl\dBi. Another important point is that the 
Neumann derivatives of -u; on S will be zero as it was the case with T. 

Since the set T was parametrized by {{—t, f{t)) : <t < 1}) the set 
E will be given by 



If we now calculate the dp = sup^g(o,p) ^^^^ that 



l + ({)2-l 1/1 



t 



sup I — sup — I ~ o'^P^- 



On the other hand 



Sp := supbJw - \y\ = Sp2 
then we need to prove that 



hminfp_^o ^f' < +oo. 
P'^Sp 



Suppose 



This means that 



Vp{x,y) = >p^Q 0. 

P(^P 

On the other hand we can find a sequence of points G S such that 
2;^ — >■ and the angle between the normal vector i^E(^fc) = (^'s(^fc)) '^K-^fc)) 
and the x-axis is large enough, i.e. 

Thus 



which is in contradiction with the fact of V\zy.\ — >■ 0. 
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Lemma 1. [Miintz Theorem] The space spanned hy {l^r""^ ^r"""^ , ...,r°'K_ 
is dense in C([0, 1]) the uniform norm if and only if 



OO 



OO. 



For a proof see [C] . 

Lemma 2. Let S = {(t,0); t E [-1,0]} and w G \^) be a 

function even in y that solves 

Aw = m Bi \ S 

(17) ij = 2t^.{-^9'o + ^) onJ: 
S = on S. 

Moreover we assume that g G C^([— 1,0]) and that g{0) = g'{0) = 0. 
Then there exist aj G M such that in polar coordinates (r, (p) 

OO 

(18) w = gq + ajr"j cos(aj<^) 

where aj are the positive solutions to 

/ \ iTc a 

(19) tan(7raj = 



Proof: From the intermediate value theorem it foUows that (19) has 
a solution aj in each interval + 1/2) (see picture |2]). Therefore it 
follows that 



■ = OO. 



CXJ - 

E- 

Also, 

^ a,- - 1 
Using Lemma [T] we can write 



OO. 



aj-l 



fi'o = «0 + E "j'^J'^ 

for some sequence of real numbers aj. By assumption gQ{0) = from 
which it follows that Sq = 0. If we define 



u = y ajr ^ cos[aj( 



simple verification to see that m is a solution to (17). 
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Next we consider v = w — u. Then v solves 

Av = in 5i \ S 
1^ = on S 

dy 

1^ = 0. 

ox 

From well known unique continuation properties of harmonic functions 
(see |AEj for instance) it follows that v =constant. The Lemma follows. 

□ 

References 

[AE] Adolfsson, Vilhelm and Escauriaza, Luis C^'" domains and unique contin- 
uation at the boundary Comm. Pure Appl. Math. Vol. 50 (1997), Nr. 10, 
935-969 

[AFP] Ambrosio, Luigi; Fusco, Nicola; Pallara, Diego Functions of bounded varia- 
tion and free discontinuity problems Oxford Mathematical Monographs. The 
Clarendon Press, Oxford University Press, New York, 2000 

[BD] Bonnet, Alexis; David, Guy Cracktip is a global Mumford-Shah minimizer 
Asterisque No. 274 (2001), vi+259 pp. 

[C] Cheney, E. W. Introduction to approximation theory McGraw-Hill Book Co. 
New York 1966 

[D] David, Guy Singular sets of minimizers for the Mumford-Shah functional 
Progress in Mathematics, 233. Birkhauser Verlag, Basel, 2005 

[KLM] Koch, Herbert; Leoni, Giovanni; Morini, Massimiliano On optimal regu- 
larity of free boundary problems and a conjecture of De Giorgi Comm. Pure 
Appl. Math. 58 (2005), no. 8, 1051-1076 

John Andersson, Mathematics Institute, University of Warwick, 
Coventry, CV4 7AL, United Kingdom 

E-mail address: J . E . AnderssonOwarwick .ac.uk 

Hayk Mikayelyan, Max-Planck-Institut fur Mathematik in den Natur- 

WISSENSCHAFTEN, InSELSTRASSE 22, 04103 LEIPZIG, GERMANY 

E-mail address: haykQmis.mpg.de 



